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Abstract 

The effect of the heavy 6-quark mass on the two, three and four-jet rates is 
studied using Lep data collected by the Delphi experiment at the Z peak 
in 1994 and 1995. The rates of 6-quark jets and light quark jets [i = uds) 
in events with n = 2, 3, and 4 jets, together with the ratio of two and four- 
jet rates of 6-quarks with respect to light-quarks, i?^^, have been measured 
with a double-tag technique using the Cambridge jet-clustering algorithm. A 
comparison between experimental results and theory (matrix element or Monte 
Carlo event generators such as Pythia, Herwig and Ariadne) is done after 
the hadronisation phase. 

Using the four-jet observable Rf, a measurement of the 6-quark mass using 
massive leading-order calculations gives: 

mb{Mz) = 3.76 ± 0.32 (stat) ± 0.17 (syst) ± 0.22 (had) ± 0.90 (theo) GeV/c^ . 

This result is compatible with previous three-jet determinations at the Mz 
energy scale and with low energy mass measurements evolved to the Mz scale 
using QCD Renormalisation Group Equations. 



(Accepted by Eur. Phys. J. C ) 



ii 



J.Abdallah^e, P.Abreu^^, W.Adam^s^ P.Adzic^^, T.Albrecht^*, R.Alemany-Fernandez^, T.AUmendingeri*, P.P.AUport^* 
U-Amaldi*", N.Amapane**, S.Amato^^, E.Anashkin^''', A.Andreazza^^, S.Andringa^^, N.Anjos^^, P.Antilogus^^ 
W-D.Apcll**, Y.Arnoudl-'^, S.Ask'', B.Asman*'', J.E.Augustin^^, A.Augustinus^, P.Baillon'', A.Ballcstrcro'*'' 
P.Bambade^i, R.Barbier^*, D.Bardini'^, G.J.Barker^^, A.Baroncelli'*", M.Battaglia^, M.Baubillier^"^, K-H.Becks^* 
M.Begalli^, A.Behrmann^*, E.Ben-Haim^i, N.Benekos^^, A.Benvenuti^, C.Berat^^, M.Berggren^^, D.Bertrand^ 
M.Besancon"! , N.Bcsson*!, D.Blochi", M.Blom^^, M.Bluj^e, M.Boncsini^", NLBoonckamp"!! , P.S.L.Bootht^* 
G.Borisov^^, O.Botner^s, B.Bouquet^^ T.J.V.Bowcock^", I.Boyko^'', M.Bracko'*'*, R.Brenner^s, E.Brodet^'' 
P.Bruckmani^, J.M.Brunet®, B.Buschbeck^^, P.Buschmann^*, M.Calvi^", T.Camporesi^, V.Canale^^, F.Carena^ 
N.Castro^S, F.CavalloS, M.Chapkin-*^, Ph.Charpcnticr", P.Chccchia''^'', R.ChiericiS, P.Chliapnikov'13, J.Chudoba'' 
S.U.ChungS, K.Cieslikis, P.CoUinsS, R.Contrii", G.Cosme^i, F.CossuttiSO, M.J.CostaS*, D.Crennell^*, J.Cuevas^s 
J.D'Hondt^, T.daSilva52, W.Da Silva^s, G.Delia Ricca^o, A.De Angelis^i , W.De Boeri*, C.De Clercq^, B.De Lotto^i 
N.De Maria**, A.De Min^'', L.dc Paula^^, L.Di Ciaccio^S, A.Di Simonc^O, K.Doroba^'^, J.Drccs^S.S, G.Eigen'' 
T.Ekelof^^, M.EUert'^s, M.Elsing^, M.C.Espirito Santo^^, G.Fanourakis^^, D.Fassouliotis^^'^, MJeindt^*^, J.Fernandez*^ 
A.Ferrer^*, F.Ferro", U.Flagmeyer^s, H.Foeth^, E.Fokitis^^, F.Fulda-Quenzer^i , J.Fuster^*, M.Gandelman^^ 

C. Garcia54, ph.GaviUet^, E.Gazis^^, R.Gokicli^'Se, B.Golob44'46, G. Gomez- Ccballos^z, P.Goncalvcs^S, E.Graziani^O 
G.Grosdidier^i, K.GrzelakS^, J.Guy^'^, C.Haagi», A.Hallgren^a, K.HamacherS*, K.Hamiltona^^ S-Haug^^*, F.Hauleris 
V.Hedberg^'', M.Hennecke", H.Herrts, J.HoffmanBe, S-O.Holmgren*'', P.J.Holt^, M.A.Houlden^*, J.N.Jackson^* 
G. Jarlskog'^'^, P.Jarry'*^, D. Jeans''^, E.K. Johansson*'^, P.Jonsson'^®, C.Joram^, L.Jungcrmann^*, F.Kapusta'^^ 
S.Katsanevas^s, E.Katsoufis^^, G.Kernel-*-*, B.P.Kersevan'**'-**', U-Kerzel^*, B.T.King^-*, N.J.Kjaer^, P.Kluit^^ 
P.Kokkinias^^, C.Kourkoumelis^, O.Kouznetsov^''', Z.Krumstein^''', M.Kucharczyk^^, J.Lamsa^, G.Leder^^, F.Ledroit^^ 
L.Lcinoncn'*'', R.Lcitncr^^, J.Lemonnc-^, V.Lcpclticrt^l ^ TXcsiak^^, W.Licbig''*, D.Liko''''', A.Lipniacka*'', J.H.Lopes^^ 
J.M.Lopez^s^ D.Loukas^^, P.Lutz"*!, LXyons^**, J-MacNaughton^s, A.Malek^*, S.Maltezos^^, F.Mandl'^s, J.Marco'*^ 
R.Marco''^, B.Marechal^^, M.Margoni^'', J-C.Marin^, C.Mariotti^, A.Markou^^, C.Martinez-Rivero*^, J.Masik" 
N.Mastroyiannopoulogl^, F.Matorras''^, CMattcuzzi^", F.Mazzucato^'', M.Mazzucato^'^, R.Mc Nulty^'*, C.Mcroni^^ 

E. Migliore"***, W.Mitaroff''^, U.Mjoemmark^^, T.Moa'*'', M.Mochi*, K.Moenig^'ii, R.Monge^*, J.Montenegro^^ 

D. Moraes^^, S.Moreno^^, P.Morettinii*, U.Mueller^s, K.Muenich^s, M.Mulders^^, L.Mundim'', W.Murray^®, B.Muryn^o 

G. Myatt^S, T.Myklebust^*, M.Nassiakou^^, F.Navarria^, K.NawrockiS", R.Nicolaidou'*! , M.Nikolenkol'''l° 
A.Oblakowska-Mucha^°, V-Obraztsov**^, A.Olshevski^'', A.Onofre^^, R.Orava^'', K.Osterberg^'^, A.Ouraou'*'^ 
A.Oyanguren^*, M.Paganoni^o, S.Paiano^, J.P.Palacios^*, H.Palkai^, Th.D.Papadopoulou^^, L.Pape^, C.Parkes^^ 

F. Parodi", U.Parzefall", A.Passeri'l°, O.Passon^S, L.Peralta^^, V.Pcrepelitsa^'l, A.Pcrrotta^, A.Pctrolinil'', J.Piedra*^ 
L.Pieri^o, F.Pierre"!!, M.Pimenta^s, E.Piotto'', T.Podobmk'14,46 ^ V.Poireau^, M.E.Pol", G.Poloki^, V.Pozdniakovi'^ 
N.Pukhaeva", A.PuUia^o, J.Ramesi^, A.Read^*, P.RebecchiS, J.Rehni*, D.Rcid^^, R.Reinhardt^s, P.Renton^s 
F.Richard^l, J.Ridky^^, M.Rivero^^, D.Rodriguez-*^, A. Romero'**, P.Ronchese^^, P.Roudeau^l, T.RovcUi^ 
V.Ruhlmann-Kleider-*!, D.Ryabtchikov'*^, A.Sadovsky*'^, L.Salmiio, J.Salt^-*, CSander**, A.Savoy-Navarro^e 
U.SchwickerathS, R.Sekulin^*, M.Siebel^®, A.Sisakiani'', G.Smadja^*, O.Smirnova^'', A.Sokolov*^^ A.Sopczak^^ 
R.Sosnowski-'''^, T.Spassov^, M.Stanitzkil* , A.Stocchi^i, J.Strauss"", B.Stugu'*, M.Szczckowski^^, M.Szcptycka''''^ 
T.Szumlak^", T.TabarcUi^", F.Tegenfeldt^^ , J.Timmermans^^, L.Tkatchevi^, M.Tobm24, S.Todorovova^^ , B.Tomc^^ 
A.Tonazzo^o, P.Tortosa^", P.Travnicek", D.Treille^, G.Tristram*, M.Trochimczuk^s, C.Troncon^s, M-L.Turlucr"*! 
I.A.Tyapkinl'', P.Tyapkin^'', S-Tzamarias^^, V.Uvarov-*^, G.Valcnti^, P.Van Dam^^, J.Van Eldik", N.van Remortcli'' 
LVanVulpen'-*, G.Vegni^s, F.Veloso^^, W.Venus^*, P.Verdier^*, V.Verzi^«, D.Vilanova^i, L.Vitale^u, V.Vrbai^ 

H. WahlenSS, A.J.Washbrook^", C.Weiseri*, D.Wicke^, J.Wickens^, G.Wilkinsoir^s, M.Wintcri°, M.Witeki'' 
CYushchenko^s, A.Zalewska^^, RZalewski^s, D.Zavrtanik^B, V.Zhuravlovi^, N.I.Zimini'^, A.Zintchenkoi'', M.Zupan^^ 



^Department of Physics and Astronomy, Iowa State University, Ames lA 50011-3160, USA 
^IIHE, ULB-VUB, Pleinlaan 2, B-1050 Brussels, Belgium 

^Physics Laboratory, University of Athens, Solonos Str. 104, GR-10680 Athens, Greece 
^Department of Physics, University of Bergen, AUcgaten 55, NO-5007 Bergen, Norway 
^Dipartimento di Fisica, Universita di Bologna and INFN, Via Irnerio 46, IT-40126 Bologna, Italy 
^Centro Brasileiro de Pesquisas Ffsicas, rua Xavier Sigaud 150, BR-22290 Rio de Janeiro, Brazil 
'^Inst. de Ffsica, Univ. Estadual do Rio de Janeiro, rua Sao Francisco Xavier 524, Rio de Janeiro, Brazil 
^College de France, Lab. de Physique Corpusculaire, IN2P3-CNRS, FR-75231 Paris Cedex 05, France 
^CERN, CH-1211 Geneva 23, Switzerland 

lOlnstitut de Rccherchcs Subatomiques, IN2P3 - CNRS/ULP - BP20, FR-67037 Strasbourg Cedex, France 
^-'Now at DESY-Zeuthen, Platanenallee 6, D-15735 Zeuthen, Germany 

i^Institute of Nuclear Physics, N.C.S.R. Dcmokritos, P.O. Box 60228, GR-15310 Athens, Greece 

i-^FZU, Inst, of Phys. of the C.A.S. High Energy Physics Division, Na Slovance 2, CZ-182 21, Praha 8, Czech Republic 
^''Dipartimento di Fisica, Universita di Genova and INFN, Via Dodecanese 33, IT-16146 Genova, Italy 
l^Institut des Sciences Nucleaires, IN2P3-CNRS, Universite de Grenoble 1, FR-38026 Grenoble Cedex, France 
^^Helsinki Institute of Physics and Department of Physical Sciences, P.O. Box 64, FIN-00014 University of Helsinki, 
Finland 

Joint Institute for Nuclear Research, Dubna, Head Post Office, P.O. Box 79, RU-101 000 Moscow, Russian Federation 
^^Institut fUr Experimentelle Kernphysik, Universitat Karlsruhe, Postfach 6980, DE-76128 Karlsruhe, Germany 
i^Institute of Nuclear Physics PAN,U1. Radzikowskiego 152, PL-31142 Krakow, Poland 

^"Faculty of Physics and Nuclear Techniques, University of Mining and Metallurgy, PL-30055 Krakow, Poland 

2iUniversite de Paris-Sud, Lab. de I'Accelerateur Lineaire, IN2P3-CNRS, Bat. 200, FR-91405 Orsay Cedex, France 

^^School of Physics and Chemistry, University of Lancaster, Lancaster LAI 4YB, UK 

23LIP, 1ST, FCUL - Av. Elias Garcia, 14-1°, PT-1000 Lisboa Codex, Portugal 

^''Department of Physics, University of Liverpool, P.O. Box 147, Liverpool L69 3BX, UK 

^^Dept. of Physics and Astronomy, Kelvin Building, University of Glasgow, Glasgow G12 8QQ, UK 

26LPNHE, IN2P3-CNRS, Univ. Paris VI ct VII, Tour 33 (RdC), 4 place Jussieu, FR-75252 Paris Cedex 05, France 

^''Department of Physics, University of Lund, Solvcgatan 14, SE-223 63 Lund, Sweden 

28Universite Claude Bernard de Lyon, IPNL, IN2P3-CNRS, FR-69622 Villeurbanne Cedex, France 

^^Dipartimento di Fisica, Universita di Milano and INFN-MILANO, Via Celoria 16, IT-20133 Milan, Italy 

30Dipartimento di Fisica, Univ. di Milano-Bicocca and INFN-MILANO, Piazza della Scienza 3, IT-20126 Milan, Italy 

31IPNP of MFF, Charles Univ., Areal MFF, V Holesovickach 2, CZ-180 00, Praha 8, Czech Republic 

32NIKHEF, Postbus 41882, NL-1009 DB Amsterdam, The Netherlands 

^^National Technical University, Physics Department, Zografou Campus, GR-15773 Athens, Greece 
^■*Physics Department, University of Oslo, Blindern, NO-0316 Oslo, Norway 
^^Dpto. Fisica, Univ. Oviedo, Avda. Calvo Sotelo s/n, ES-33007 Oviedo, Spain 
^^Department of Physics, University of Oxford, Keble Road, Oxford OXl 3RH, UK 
^'^Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, IT-35131 Padua, Italy 
3*Rutherford Appleton Laboratory, Chilton, Didcot OXll OQX, UK 

^^Dipartimento di Fisica, Universita di Roma II and INFN, Tor Vergata, IT-00173 Rome, Italy 
^''Dipartimento di Fisica, Universita di Roma III and INFN, Via della Vasca Navale 84, IT-00146 Rome, Italy 
''iDAPNIA/Servicc de Physique des Particules, CEA-Saclay, FR-91191 Gif-sur-Yvette Cedex, France 
''^Instituto de Fisica de Cantabria (CSIC-UC), Avda. los Castros s/n, ES-39006 Santander, Spain 
''^Inst. for High Energy Physics, Serpukov P.O. Box 35, Protvino, (Moscow Region), Russian Federation 
''''j. Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia 

''^Laboratory for Astroparticle Physics, University of Nova Gorica, Kostanjeviska 16a, SI-5000 Nova Gorica, Slovenia 
''^ Department of Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
''''Fysikum, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

'** Dipartimento di Fisica Spcriiueiitalc, Universita di Torino and INFN, Via P. Giuria 1, IT-10125 Turin, Italy 
*^INFN,Sezione di Torino and Dipartimento di Fisica Teorica, Universita di Torino, Via Giuria 1, IT-10125 Turin, Italy 
5°Dipartimento di Fisica, Universita di Trieste and INFN, Via A. Valerio 2, IT-34127 Trieste, Italy 
silstituto di Fisica, Universita di Udine and INFN, IT-33100 Udine, Italy 

s^Univ. Federal do Rio de Janeiro, CP. 68528 Cidade Univ., Ilha do Fundao BR-21945-970 Rio de Janeiro, Brazil 

^^Department of Radiation Sciences, University of Uppsala, P.O. Box 535, SE-751 21 Uppsala, Sweden 

s^IFIC, Valencia-CSIC, and D.F.A.M.N., U. de Valencia, Avda. Dr. Moliner 50, ES-46100 Burjassot (Valencia), Spain 

^^Institut fiir Hochenergiephysik, Osterr. Akad. d. Wissensch., Nikolsdorfergasse 18, AT-1050 Vienna, Austria 

®®Inst. Nuclear Studies and University of Warsaw, Ul. Hoza 69, PL-00681 Warsaw, Poland 

^''Now at University of Warwick, Coventry CV4 7AL, UK 

^*Fachbereich Physik, University of Wuppertal, Postfach 100 127, DE-42097 Wuppertal, Germany 
t deceased 



1 



1 Introduction 

Mass corrections to the Z^bb coupling are of order (m^/M|), which is too small to 
be measured at Lep and Slc. For some inclusive observables, like jet-rates, the effect is 
enhanced as {ml/M'^)/ycut, where ijcut is the jet resolution parameter [1]. The effect of 
the 6-quark mass in the production of three-jet event topologies at the Z peak has for 
instance already been measured at Lep and Slc [2-5]. Multi-jet topologies with &-quarks 
appear both as signal and background in searches and precision measurements at current 
and future colliders. Their study, together with that of the gluon emission from massive 
quarks, is an effective tool to probe the fundamental short-distance QCD features of the 
Standard Model and is important to test the modelling of b and light-quark jets available 
in calculations and generators. 

This study generalizes the methods described in references [2, 6] and presents the 
measurement of the normalized n-jet production partial widths for Z-decays into 6-quark 
or light quark pairs: 

Z^qq 

q = b,e (i^uds), (1) 

depending on the ycut value of the Cambridge jet-clustering algorithm [7] which is used 
here. 

The effect of the heavy 6-quark mass on jet rates is studied by measuring the double- 
ratio observable: 

The Delphi data collected during the years 1994 and 1995 at a centre-of-mass energy 
of ^/s ~ Mz have been analysed. Experimental results are compared to the hadronic final 
state simulated by the fragmentation models of Pythia 6.156 [8], Herwig 6.2 [9] and 
Ariadne 4.08 [10] and to matrix element (ME) calculations folded with a hadronisation 
correction. Therefore, the data are corrected for detector and kinematical effects, while 
ME calculations, computed at parton level, are corrected for hadronisation. 

In order to extract the 6-quark mass information from i?^^ measurements, massive ME 
calculations performed in terms of both the pole mass M5 and the running mass ^^(/i) 
are used. Jet-rate calculations are only available to 0{p?g) [11-13], therefore massive four- 
jet observables can only be described to leading-order (LO) accuracy. The 6-quark mass 
obtained from I^l using such LO calculations is compared to the three-jet results [6] 
and to mass values at threshold [14] evolved to the Mz scale using Renormalisation 
Group Equations (RGE). An approximate massless NLO correction is also tried as an 
improvement. 

The precision of 6-mass measurements from three-jet events is limited by systematic 
uncertainties (hadronisation, 6-tagging and theory). The four-jet observable li'l has a 
larger statistical error but its sensitivity to the 6-quark mass is higher because, most 
probably, the emission of two gluons is involved. The four-jet topology thus provides 
a complementary measurement in which the systematic uncertainties can be expected 
to be partly different. In this analysis, flavour jet-rates arc measured using a double- 
tag technique which measures signal and background efficiencies from data in a self- 
calibrating way, reducing the systematics and allowing for a useful cross-check of previous 
measurements [6]. 



■n=2, 3,4(11 cut 



cut J 



tot 
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2 The DELPHI detector 

Delphi was a hermetic detector located at the Lep accelerator, with a supercon- 
ducting solenoid providing a uniform magnetic field of 1.23 T parallel to the beam axis 
throughout the central tracking device volume. A detailed description of its design and 
performance is presented in [18,19]. 

In the Delphi coordinate system, the z axis is oriented along the direction of the 
electron beam. The polar angle Q is measured with respect to the z axis, is the 
azimuthal angle in the plane transverse to the z axis and R — \J + is the radial 
coordinate. 

The main tracking devices in Delphi were the silicon Vertex Detector (VD), a jet 
chamber Inner Detector (ID) and a Time Projection Chamber (TPC). They were located 
in the immediate vicinity of the interaction region to reduce the amount of material 
between the beam and the detector. At a larger distance, the tracking was completed by 
a drift chamber Outer Detector (OD) covering the barrel region (40° < 9 < 140°) and 
two sets of drift chambers, FCA and FCB, located in the endcaps. 

The VD was the detector closest to the interaction point. In 1994 and 1995 it con- 
sisted of three coaxial cylinders, the inner and outermost ones consisting of double-sided 
detectors with orthogonal strips, allowing the measurement of both Rcj) and z coordinates. 

Electron and photon identification was provided by electromagnetic calorimeters: the 
High Density Projection Chamber (HPC) in the barrel and a lead-glass calorimeter 
(FEMC) in the endcaps. Hadronic energy was measured in the hadronic calorimeter 
(HCAL). 

3 Data analysis 

First, the sample of Z hadronic decays, i.e. Z^qq events was selected. Then the 
different jet-topologies were identified using the Cambridge jet-clustering algorithm [7]^, 
and b and light-quark samples were separated using the DELPHI flavour tagging methods, 
based on properties of the long-lived heavy 5-hadrons. Experimental results were then 
corrected for detector and acceptance effects in two different ways, depending on the 
observable and topology, as explained in Section 3.3. Matrix element and event generator 
predictions were corrected for hadronisation effects from the parton to the hadron level. 
The parton level is defined as the final state of the parton shower (in Pythia and 
Herwig) or dipole cascade (in Ariadne) in the simulation, before hadronisation. These 
corrections are discussed in Section 4. 

3.1 Event selection 

Total numbers of 1 484 000 and 750 000 hadronic Z boson decays, collected at the Z 
resonance by Delphi during the years 1994 and 1995, respectively, have been analysed 
in order to study mass effects in multi-jet topologies^. 

Hadronic events were selected in the same way as in reference [2] (see Table 1, left): 

• Charged and neutral particles were reconstructed as tracks and energy depositions 
in the detector. A first selection was applied to ensure a reliable determination of 
their momenta and energies; 

^In our analysis, the values of the ordering and resolution parameters were taken to be equal, Vij = yij . 
^Earlier data samples were not considered as the VD setup was less complete and resulted in a less precise flavour 
identification, which is crucial for this analysis. 
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• The information from the accepted tracks was combined event-by-event and hadronic 
events were selected according to global event properties. 

Finally, a total sample of 1 150 000 Z hadronic decays was selected. Then jets were 
reconstructed with the Cambridge algorithm. In order to reduce the impact of particle 
losses and wrong energy-momentum assignment to jets, further kinematical selections 
were applied, which were slightly different for each jet topology (see Table 1, right). 

Simulated events were produced with the DELPHI simulation program Delsim [19], 
based on Pythia 7.3 tuned to Delphi data [20], and were then passed through the same 
reconstruction and analysis chain as the experimental data. The simulated events were 
reweighted in order to reproduce the measured rates of \h and cc-quark pairs arising from 
the gluon splitting processes [21] {g^-^ = 0.00254 ± 0.00051, gcc = 0.0296 ± 0.0038), which 
are significantly larger than those in the standard simulation. 



Charged 

Particle 

Selection 


Pch > 0.1 GeV/c 
25° < ^ < 155° 
L > 50 cm 

d < 5 cm in R(f) plane 
d < 10 cm in z direction 


Neutral 
Cluster 
Selection 


> 0.5 GeV,40° <9< 140° 
^FEMC > 0.5 GeV,8° < ^ < 36° 
£.FEMC > 0.5 GeV,144° < ^ < 172° 
E^^^ > 1 GeV, 10° < ^ < 170° 


Event 
Selection 


Nch > 5 

Ech > 15 GeV 

IE/9i| <6,i= l:...,iV,, 

No particle with p,,i^ > 10 Gc\Yc 



2-jct 


45° < Othrust < 135° 


3-jet 


45° < Othrust < 135° 
Nf' > 1 
Ej > 1 GeV 
25° < 9j < 155° 
> 359°, 2<j 


4-jet 


32° < Othrust ^ 148° 

Ej > 1 GeV, 
25° < Oj < 155° 



Table 1: (Left) Particle and event selections: Pch is the momentum of charged particles, 
L their measured track length, d their impact parameter with respect to the interaction 
point and their charge, E^i is the energy of neutral clusters in the calorimeters, N^h 
is the number of charged particles and E^h their total energy in the event. (Right) 
Kinematical selections for jets in accepted events: Othrust is the polar angle of the thrust 
of the event, Nj'^ the charged multiplicity in the jet, Ej the jet energy and Oj the angle 
between the jet and the beam axis. For three-jet events, an additional planarity cut is 
applied on the sum of all jet pair angles, (pij. 



3.2 6-tagging 

The identification of 6-quark events in Delphi was based on the properties of a B- 
hadron such as its large mass and the large impact parameter of its decay products. A 
jet estimator variable Xj^t was built as an optimal combination of five discriminating 
variables [22]. The most discriminant one was the probability of having all charged 
particles in the jet produced at the event interaction point. The use of this variable 
alone defined the impact-parameter technique. The additional variables were used only 
when a secondary vertex (SV) was reconstructed. These variables were, for all particles 
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attached to the SV: the invariant mass, the fraction of the charged jet energy, the sum 
of all transverse momenta and the rapidity of each particle. The information from all 
five variables was combined into a single estimator Xj^t in an almost optimal way which 
provided discrimination between heavy and light jets with high purity and efficiency. To 
obtain 6(light)-quark enriched samples, jets with an estimator value above (below) a given 
threshold Xj^t > X^^^ {^jet < ^jet) were selected. To tag events, the value of the two 
highest 6-tagging jet variables were combined into an event estimator, X^^ = Xj^^ + X^^^. 

In the present analysis, this approach has been associated to a double-tag tech- 
nique [23], which measures flavour-tagging efficiencies directly from data. 

Using the two jets with highest 6-tagging variables as the flavour jets (jets which 
are expected to contain a primary quark) makes no distinction between primary quarks 
originating in the Z decay and secondary production of b and c-quarks from gluons 
(g^bb,cc), a process referred to as gluon splitting and which constitutes a significant 
part of the systematic uncertainty in multi-jet flavour-observables (sec Section 4.3). 

To reduce the sample contamination from gluon splitting in four-jet events, the flavour 
jets were defined as follows: the most energetic jet in each event is identified as the first 
flavour jet. Remaining jets are ordered by angular proximity to it. The closest jet is 
discarded making the hypothesis that it is a gluon coming from the same primary quark. 
The second 6(light)-flavour jet is that with the highest 6-tag (lowest 6-tag) estimator 
among the two remaining jets. In this way, energy and angle information is combined to 
define the fiavour-jets. As an additional selection, an event is not classified as bb if the 
most 6-tagged jet is not among the two most energetic jets; this last selection reduces the 
uncertainty from (7^,^ and Qcc by a factor two. The effect of the remaining contamination 
due to gluon splitting is included in the gluon-splitting uncertainty and is well below the 
statistical uncertainties (Tables 4 and 5). 

3.3 Overview of the correction method 
3.3.1 Event-tag 

To correct the two-jet observable i?2^ for detector effects and the flavour tagging pro- 
cedure, the event-tag method described in reference [2] was used: 



where the measured rate R2 is corrected by using purities of the inclusive samples, 
Cq = Nq/Nq (the fraction of qq events tagged in the Q category), and detector corrections 
taken from the Delsim simulation, = R^q/^I (where R^q is the two-jet rate of qq 
events tagged as Q). The factor Rf = R^R^ is taken from the simulation. Table 2 
summarizes the number of events selected in the 1994 data in each flavour sample for the 
chosen working points of purity = = 98% {Pl = c^l — 73%, L = uds) and efficiency 
of = 38% (e^ = 58%) for 6-flavoured (light-flavoured) events (where = NQ/Ng, the 
ratio of tagged events of a given flavour to the total number of events of the same flavour), 
respectively. 

The event-tag method has the advantage of applying the flavour-tagging procedure 
only in the inclusive sample, before events are classified into jet topologies. 



„b TDM-det _ h Jb 
<-L"2L-'T'2 <-b"2B 



M-det 




(3) 
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3.3.2 Double-jet tag 

The event-tag method, if the jet sample is topologically very different from the inchisive 
one, can introduce important biases. To prevent this, in the B}f measurement 6-tagging 
is applied to jets. The observable in Eq. 2 is rewritten as: 

^ ( T{Z^a) \ Nl/N, ^ n-R,-R \ Nl 
' \T{Z^hh)) Ni/N, V Ri> ) Ni' ^ > 

The global normalisation can be obtained directly from the world average values of and 
Rc [14]: 

Rb = 0.21629 ±0.00066, (5) 
Rc = 0.1721 ± 0.0030, 

which implies a ±6%o uncertainty on Ft^. A double-tag technique is used: the total 
number of four-jet events, A^4, the corresponding numbers for a given flavour Nl, q — 

b, udsc, and the tagging efficiencies and ^if^sc ^'^^ obtained from comparing the number 
of four-jet events where two jets are tagged as h or udsc to the number of events where 
a single jet is tagged. This is done by solving the following set of equations: 

M = {Nlel + {N,-Nl)er-'], (6) 

\m,b = {Nlele'B + {N,-Nl)er-'eT-']. (7) 

U,BB = {Nle\e'^B + {N,-Nt)er-'el^-'], (8) 

and equivalent equations for the udsc-tagged samples. The left hand side of these equa- 
tions are the measured quantities. A4 is the number of measured four-jet events. For 
each event the two jets which are most likely to contain a primary quark (flavour jets, see 
above) are selected and the flavour identiflcation is done independently for both jets: Mae 
is the number of jets tagged as B (with a maximum possible value of 2A4, two from each 
event) and A^ss is the number of events where the two flavour jets are simultaneously 
tagged as B. With this method, the jet-rates R\ and i?4 are measured independently, 
together with the efficiencies and ^^osc- To accomplish this, double-jet tagging effi- 
ciencies Eqq are related to the single jet-tagging efficiencies through correlation factors 
defined from egg = N^q/N'^ = egeg(l -|- pg). Here, charm-events have been included 
in the udsc-tagged category: the hght-quark content is extracted from TV"'^*'^ after 
dividing by a factor (1 + N^/N^) obtained from Monte Carlo event generators. Only 
hadronic event-selection efficiencies for each flavour, = N^'^'' /Ng, mistagging efficien- 
cies, eg"*^"^ and Cgg "^, and flavour correlations for b and light-tagging are computed from 
the simulation. 

This procedure can be easily generalised to cover n — 2, 3-jet topologies in order to 
measure both jct-ratcs (i?^, R^) and the double-ratios (R^^^) independently. Due to the 
6%o uncertainty from the global normalisation, the double-tag measurements for i?2^3 
are less precise than the corresponding event-tag result. However, they serve as a useful 
cross-check both of the final result and on the consistency between data and simulation 
for the flavour-tagging efficiencies. 

Results with this method have a better stability with respect to the value of the flavour- 
tagging threshold, and are more consistent with each other"^. The flavour composition of 
the 1994 sample is shown as an example in Tables 2 and 3. The stability obtained in the 
case of the four-jet rates is shown in Figure 1 for the 1994 and 1995 data samples. 

^ From the relation = + N^dsc^ tj^e double-ratio can be obtained independently in two ways, starting either 
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Event-tag method 


Flavour 


Inclusive 


2 jets 


cut 


Purity 


Efficiency 


B 


111440 


75147 


Xe. > 1.10 


98% 


38% 


L 


678282 


414912 


Xev < 0.40 


73% 


58% 



Table 2: Flavour composition of the 1994 sample {ycut = 0.0065). The number of events in 
the inclusive and 2-jet samples are shown separately for B and L — uds tagged events for 
the chosen flavour-tagging working points. Purity and efficiency are also shown. Similar 
numbers were found with the 1995 data. 



double-tag method 



Topology 


Q 


QQ 


cut 


Purity 


Efficiency 


2 jets [B] 
2 jets (I) 


136228 
610757 


35187 
213716 


X^r, > +0.33 
X2, < -0.92 


92% 
95% 


57% 
78% 


3 jets {B) 
3 jets (L) 


66034 
362396 


15356 
147605 


Xsj > +0.19 
X:^.j < -0.64 


87% 
93% 


53% 
84% 


4 jets {B) 
4 jets (L) 


10720 
91042 


2191 
36773 


Xij > +0.05 
X^j < -0.64 


84% 
89% 


35% 
86% 



Table 3: Flavour composition of the 1994 sample 0.0065) tagged as n-jet 6-quark 

(B) and udsc-quark events (L) for the different jet topologies analysed, n = 2, 3 and 
4 jets. Four-jet tagging uses the method described in Section 3.2 for the definition of 
flavour jets. Similar numbers were found with the 1995 data. 

4 Results 

The single-flavour jet rates i?^, n = 2,3,4-jets, and the four-jet observable Rf, are 
measured with the double-tag technique, while the two-jet observable i?2^ is measured us- 
ing the event-tag method described in [2] . A description of the experimental uncertainties 
considered in the analysis is given in Section 4.3. Theoretical uncertainties, arising in 
the comparison between ME predictions and the four-jet observable, are discussed in 
Sections 4.5 and 4.6. 

4.1 Single jet-rates, 

The measured rates {n = 2,3,4 jets, q = b ot £ = uds) are shown in Figure 2a 
together with predictions from the Pythia 6.156, Herwig 6.2 and Ariadne 4.08 gen- 
erators tuned to Delphi data [20] (see Section 4.6 for the choice of the 6-quark mass 
parameter in the generators). The detailed breakdown of the uncertainties of the mea- 
sured jet-rates is shown in Table 4. The measurements in 1994 and 1995 were found 
to be incompatible with each other at the two standard deviations level, indicating that 
some systematic effect was not taken into account in the three-jet light-quark rate. The 
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systematic tagging uncertainty in R\ was increased in order to fully cover this difference. 
Only the uncertainty in i?3 was increased since the 6-tagging was developed from 2-jet 
events yielding reliable i?2 results, and in 4-jet events the 6-tagging applies different cuts 
on angle and energy. The consistency of the experimental results and the prediction from 
the three event generators is shown in Figures 2b-c: the Herwig 6.2 and Ariadne 
4.08 generators provide a reasonable description of the six observables in the region of 
Ucut between 0.001 and 0.010. Pythia 6.156 gives the best description of R\, but is 
inconsistent with the other jet-measurements at the three standard deviations level. 





i?2 


i?2 


i?3 




i?4 




Value 


0.6224 


0.6034 


0.3004 


0.3150 


0.0598 


0.0676 


Statistical (data) 


±0.0019 


±0.0008 


±0.0016 


±0.0006 


±0.0007 


±0.0004 


Statistical (sim.) 


±0.0012 


±0.0005 


±0.0009 


±0.0004 


±0.0006 


±0.0003 


Tagging 


±0.0004 


±0.0008 


±0.0006 


±0.0025 


±0.0001 


±0.0002 


Normalisation 


±0.0018 


±0.0030 


±0.0009 


±0.0016 


±0.0002 


±0.0003 


9bb 


±0.0003 


< 0.0001 


±0.0009 


< 0.0001 


±0.0003 


< 0.0001 


9cc 


±0.0006 


±0.0002 


±0.0010 


< 0.0001 


±0.0002 


< 0.0001 


Total systematics 


±0.0020 


±0.0031 


±0.0017 


±0.0030 


±0.0004 


±0.0004 


Total statistical 


±0.0023 


±0.0009 


±0.0018 


±0.0007 


±0.0009 


±0.0005 


Total uncertainty 


±0.0030 


±0.0033 


±0.0025 


±0.0031 


±0.0010 


±0.0006 



Table 4: Breakdown of uncertainties for the i?^ jet-rate measurements at a reference 
Ucut — 0.0065. The definition of each systematic contribution is given in Section 4.3. 



4.2 Double-ratios, i?^^ 

The measured double-ratios (n = 2,3,4 jets) are shown in Figures 3-4 together 
with predictions from the Pythia 6.156, Herwig 6.2 and Ariadne 4.08 generators 
tuned to Delphi data [20] (see Section 4.6 for the choice of the 6-quark mass parameter 
in the generators). 

Results for R^^ and from the event-tag and double-tag methods are shown in 
Figure 3 (event-tag results for Rf are taken from [6]). i?2^ is not described well by 
either of the generators in the full Ucut range. In all cases, both methods give consistent 
results within one standard deviation. A better experimental precision is found with 
the event-tag, because the global normalisation uncertainty is absent in this case and 
because flavour-tagging uncertainties cancel to first order in the products Cgrf^g (see 
Eq. 3). Statistical uncertainties in the event-tag result are also smaller, as more data 
events are considered and as statistical fluctuations are partially reduced in the ratios 
of the jet and inclusive samples. The detailed breakdown of the uncertainties of the 
measured double-ratios is shown in Table 5 for the event-tag method. 

The Rf result with the double-tag method is shown in Figure 4a, while the experi- 
mental systematics breakdown is summarized in Table 5. At ycut values above 0.004 the 
measurement is dominated by statistical uncertainties, while for very low values of ycut the 
data samples increase and the global normalisation uncertainty dominates. Gluon split- 
ting uncertainties are kept low in the whole ycut range thanks to the dedicated anti-gluon 
splitting cut (see Section 3.2). Herwig provides the best description, being compatible 
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LO (GeV/c") 


NLO (GeV/c") 




K2 


K [6] 




mb{Mz) = Mb 


mb{Mz) 


Mb 


Value 


1.0440 


0.9570 


0.883 


3.76 


3.46 


5.07 


Statistical (data) 


±0.0021 




±0.012 


±0.25 


±0.27 


±0.35 


Statistical (sim.) 


±0.0012 




±0.010 


±0.20 


±0.22 


±0.28 


Tagging 


±0.0009 




±0.003 


±0.07 


±0.08 


±0.10 


Normalisation 


itU.UUUO 




itU.UUO 


1 r\ 1 1 
±0.11 




1 n 1 p. 


9bb 


±0.0007 




±0.005 


-i-D 1 n 


±0.10 


±0.13 


Qcc 


±0.0003 




±0.003 


±0.06 


±0.06 


±0.08 


Total systematics 


±0.0013 


±0.0027 


±0.008 


±0.17 


±0.19 


±0.24 


Total statistic-al 


±0.002 1 


±0.0037 


±0.015 


±0.32 


±0.35 


±0.16 


T( )tal exi)eriiiieiital 


±0.0027 


±0.0016 


±0.017 


±0.36 


±0.10 


±0.52 


Modelling 








±0.22 


±0.24 


±0.32 


Theoretical 








±0.90 


±0.44 


±0.57 


Total uncertainty 


±0.0027 


±0.0046 


±0.017 


±0.99 


±0.64 


±0.83 



Table 5: Breakdown of uncertainties for the i?^^ (n = 2,3,4) double-ratio measure- 
ments. The three-jet result is taken from [6] and shown here for completeness. The 
two and three-jet measurements are based on the event-tag method, while R!f uses the 
double-tag technique as explained in Section 3.3. The 6-mass values (running and pole) 
extracted from R^^ at reference ycut = 0.0065 are also shown, both for the massive LO 
{Mf, = mb{Mz) ) and approximate NLO calculations. Experimental and modelling un- 
certainties (experimental tuning and hadronisation model in the simulation) are detailed 
separately. 

with the experimental data in the whole ycut range. However, the Pythia prediction 
is only 1.5 standard deviations away in the large ycut region; Ariadne provides a good 
description of the data in the region ycut > 0.005, while for lower values of ycut it tends 
to underestimate the mass effect. 

4.3 Experimental uncertainties 

Experimental uncertainties arise in the process of correcting the detector-level mea- 
surement to hadron level, and are due to imperfections in the physics and detector mod- 
elling in the Delsim simulation used in the correction procedure. The following sources 
have been considered in this analysis: 

• Statistical: these uncertainties are due to the limited size of the experimental and 
simulated data samples. They are estimated from a toy simulation based on Poisson 
statistics. Central values were taken from the data and simulated samples, and 
correlations between the different quantities were accounted for by building up the 
corresponding covariance matrix. 

• Gluon splitting: the identification of primary 6-quarks is based on the presence of 
long-lived B and D-hadrons in the final state. However, light-quark events with 
gluon radiation splitting into secondary heavy quarks can produce a similar signa- 
ture. The correction procedure is very sensitive to the gluon splitting rates in the 
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Monte Carlo simulation through the signal and background efficiencies [6]. Their 
value was varied in the range of their quoted uncertainties [21] and the observed 
change in the observables was added in quadrature and taken to represent the cor- 
responding uncertainty. 

• Normalisation: the uncertainty on the global normalisation Rb/Re is estimated by 
varying the world average values of Rh and Ri = (1 — i?^ — Re) in the range of their 
quoted uncertainties [14], and taking the maximum variation in the final observable 
as the global normalisation uncertainty. This results in a 6%o relative uncertainty 
and is y^^ut independent. The uncertainty from the charm-/light-quark normalisation 
factor {R^ = R^/ R^^) is estimated as half the maximum difference obtained by using 
as input to the measurement the prediction from the three event generators used: 
Pythia 6.156, Herwig 6.2 and Ariadne 4.08. 

• Flavour-tagging: signal efficiencies (e^^ and eJJ^*'^) arc measured from data and there- 
fore do not contribute to the total uncertainty for the double-tag technique. To 
estimate the uncertainty due to the imperfect description of background efficiencies 
and flavour correlations (p^g) in the simulation, the calibration of the 6-tagging 
in the simulation was exchanged with the calibration obtained from data, which 
gives a poorer description of the lifetime probability [23]. Twice the observed differ- 
ence was conservatively taken as the flavour-tagging uncertainty. For the event-tag 
technique, the related uncertainty was estimated as in [6] by varying the tagging ef- 
ficiencies within their uncertainties: Ae^g/e^g = 3% and Ae^^/e^^ = 8% evaluated 
in reference [23]. The effect of mistagging efficiency was estimated by considering 
light-tagging as equivalent to anti 6-tagging, i.e. Ae^^ = Ae^^ for q = b,c,£ for the 
same cut value. 

4.4 Hadronisation corrections 

To compare parton-level fixed order ME calculations of i?4^~^"'^* with experimental 
results, they must be corrected for hadronisation effects: 

< = HfR^t""''. (9) 

The corrections H\^{ycut) relating parton to hadron observables are taken to be linear 
bin-to-bin factors. 

Three different generators, each tuned independently to the Delphi data [20], were 
used in this analysis: Pythia 6.156, Herwig 6.2 and Ariadne 4.08. It was found 
that the Herwig and Ariadne event generators are consistent both with the theoreti- 
cal predictions at the parton level (within the theoretical uncertainty) and the data (see 
Figure 4) for a large range of Ucut- The hadronisation corrections computed with the 
three generators are shown in Figure 4b. The average of the Herwig and Ariadne pre- 
dictions was used to correct the massive ME theoretical calculations (in the region of 
Ucut studied here, the hadronisation correction computed from Pythia is contained in 
the band defined by the Herwig and Ariadne corrections). 

4.5 6-quark mass extraction and approximate NLO ME calcu- 
lation 

For a given flavour the n-jet rate is dcflned as the normalised n-jet cross-section 
-^n — \^n/'^tot\'^^'*^ ■ Theoretically, it is convenient to use the double-ratios = R\J R^^ 
as in this observable most of the higher-order electroweak corrections, the first order 
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dependence on and, to some extent also neglected higher-order terms in a^, cancel 
out. Massive ME theoretical calculations exist up to order [11-13] and describe the 
2, 3 and 4-jet rates for heavy (6, c) and light quarks {i — uds). Such calculations, when 
performed in the on — shell scheme in terms of the pole mass, Mq, can be rewritten in 
terms of the running mass, ruq, defined in the MS scheme, using the following order 
as relation: 



Both mass definitions are equivalent at LO (see Eq. 10). For ycut — 0.0065, a value 
within a region with good stability, high sensitivity and small hadronisation corrections, 
the following 6-quark mass value was obtained: 

Mf, = mb{Mz) = 3.76 ± 0.32 (stat) ± 0.17 (syst) ± 0.22 (had) ± 0.90 (theo) GeV/c^. 

The theoretical uncertainty is estimated as half the difference between the LO pre- 
diction for the running and pole 6-quark mass definitions (see Figure 4b). 

To extract a meaningful 6-quark running mass from the four-jet observable by means 
of Eq. 10, the NLO correction to i?^ would be needed, which is only available for massless 
quarks [15, 16]. However, an improvement of the LO estimation can be obtained if most 
of the mass effect is contained in the LO term and hence the NLO correction to Rf can 
be approximated as massless [17]: 



where the LO functions A^,A^ are taken from [11-13] and the NLO massless term 
from [15, 16]. As for the case of Rf [24], it was found that: 

• the NLO corrections using the pole and running mass definitions were both within 
the uncertainty band defined by the two LO curves; 

• the running mass definition results in a smaller correction at NLO than the pole 
mass. 

The 6-mass values obtained from R^ using this approximation are shown in Figure 5b-c. 
They are found to be stable in the region i/cut > 0.003 and consistent with mass results 
obtained from R'^^ (both at LO and NLO) and predicted values from QCD calculations at 
low energy evolved to Mz using the RGE. For the running mass calculation, the massless 
NLO correction is small and results in very little effect. On the contrary, for the pole 
mass the NLO correction is about 10%, leading to sizeable effects. 

For the running 6-quark mass definition, the theoretical prediction of R'f is taken to 
be the central value of the following, in principle equivalent, four calculations: (a) Full 
ratio as in Eq. 11, expressed in terms of the running mass by means of Eq. 10 at the scale 
jj, — Mz; (b) Same, but using Eq. 10 at an arbitrary scale //q = Mi, and evolving the result 
to /i = Mz via the RGE to obtain mi,{Mz)', (c) Series expansion of Eq. 11, expressed 
in terms of mb{Mz) as in the first method; (d) Same, but introducing an arbitrary 
intermediate scale as in the second method. The pole mass prediction is obtained in a 
similar way. The resulting predictions for R^ are shown in Figure 4b for a reference 
6-quark mass obtained by evolving the average of low energy measurements 7715(7725) = 
4.20 ±0.07 GeV/c^ [14] to the Mz scale, mb{Mz) = 2.84 ±0.06 GcV/c^, or by translating 
it to a pole mass value: M^ = 4.94 ± 0.08GeV/c^. The strong coupling constant value 
used was as{Mz) = 0.1202 ± 0.0050 [25]. 




(10) 




A\mb)al + B^al 



(11) 
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4.6 Theoretical and modelling uncertainties 

The following sources of systematic uncertainty have been considered for the compar- 
ison of the corrected four-jet ME calculations with the experimental results: 

• Theoretical uncertainties, due to missing higher orders in matrix element calculations 
and to the use of massless next-to-leading corrections for the mass extraction, cannot 
be rigorously estimated in the case of four-jets. However, following a comparison 
between the same approximation applied to R^^ with the full massive calculation 
available in this case, this uncertainty was conservatively taken to be twice the max- 
imum difference between the four predictions defined in Section 4.5. The theoretical 
uncertainty is responsible for about 0.4 — 0.5 GeV/c^ in the uncertainty of the final 
result, and it is almost independent of i/cut- Although lower than in the case of the 
LO calculation, it is three times higher than in the completely massive three-jet 
calculation. 

• Modelling uncertainties, related to the correction for hadronisation effects of the 
theoretical calculations at parton level using Monte Carlo event generators. This 
includes the uncertainty on the tuned values of the free parameters in each model 
(including the 6- mass parameter entering in the parton shower [6]) and the modelling 
of hadronisation. The size of the modelling uncertainty is estimated as half the differ- 
ence between the predictions from Herwig 6.2 and Ariadne 4.08^. To include the 
6-mass uncertainty in the estimation of the hadronisation systematic uncertainty, the 
mass parameter in Herwig and Ariadne was varied within ±0.125 GcV/c^ around 
their central values in order to maximize the difference between both predictions. 
This was achieved by setting the mass parameter to = 4.85 GeV/c^ in both 
generators. The total modelling uncertainty amounts to ±(1 — 2)% in the region of 
Ucut > 0.004, corresponding to about ±0.2 GeV/c^ in terms of both the running and 
pole mass results. The contribution from varying the mass parameter amounts to 
about ±0.1 GcV/c^. 

The breakdown of the theoretical and modelling uncertainties in the 6-quark mass results 
obtained from is detailed in Table 5. 

5 Summary and conclusions 

A new determination of the hadron-level i?^ and R^^"^^ jet-rates {n = 2, 3, 4 jets) has 
been performed, using flavour tagging only in each n-jet sample and obtaining the global 
normalisation of the observables from the world average Ri, and Rc measurements [14]. 
This measurement is based on a double-tag technique which measures the flavour-tagging 
efficiencies directly from data, thereby reducing systematic uncertainties. 

Double-ratio observables are also studied: R\^ is obtained from the four-jet rates R\ 
and i?4 using this double-tag technique, and i?2^ using the event-tag method defined in 
reference [6]. Results from i?2^ (and from the previous measurements of R^^ in [6]) are 
also cross-checked. 

Results are presented at hadron level, in order to allow for future comparisons without 
having to unfold hadronisation and detector corrections applied to the data (a summary 
of jet-rate results as a function of ycut is shown in Tables 6-7). They are compared to 
three Monte Carlo event generators: Pythia 6.156, Herwig 6.2 and Ariadne 4.08, 
tuned to Delphi data [20]. The Herwig 6.2 generator gives the best overall description 

"^Thc result obtained with the Pythia 6.156 event-generator is compatible with the quoted results within the modelling 
uncertainty. 
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of flavour jet-rates, i?J^ and i?^, but Ariadne 4.08 provides the best results for i?^. 
For double-ratios, Herwig 6.2 gives also the best description. However, the two-jet 
observable is not satisfactorily described by any of the three generators considered. 

A new determination of the 6-quark mass in the four-jet topology has been performed 
using the Cambridge jet-clustering algorithm [7]. The mass is measured by comparing 
the experimental results of Bfjl at Ucut = 0.0065 with fixed order ME massive LO calcula- 
tions assuming the universality of the strong coupling constant, ag- The measured value 
is: 

mb{Mz) = 3.76 ± 0.32 (stat) ± 0.17 (syst) ± 0.22 (had) ± 0.90 (theo) GeV/c^ 

A procedure to approximate the NLO corrections with the massless component in order 
to improve the result has been tested successfully with the three-jet massive calculations. 
The measured value of the running &-quark mass when applying this method to the 
four-jet observable is: 

mi,{Mz) = 3.46 ± 0.35 (stat) ± 0.19 (syst) ± 0.24 (mod) ± 0.44 (theo) GeV/c^ 

and the corresponding value for the pole mass is: 

Mb = 5.07 ± 0.46 (stat) ± 0.24 (syst) ± 0.32 (mod) ± 0.57 (theo) GeV/c^. 

These results agree within the uncertainties with the values obtained evolving the av- 
erage of low energy measurements mb{mb) — 4.20 ± 0.07 GeV/c^ [14] to the Mz scale 
using the RGE: mb{Mz) = 2.84 ± 0.06 GeV/c^, or by translating it to a pole mass 
value: Mb = 4.94 ± 0.08GeV/c^. The values of mb{Mz) obtained from the LO and 
approximate NLO R'f calculations are shown in Figure 6 together with results from 
other measurements at the Mz scale, in particular the most precise result from Rf, 
mb{Mz) = 2.85 ± 0.32 GeV/c^ [6], as well as results at low energy from semilcptonic B- 
decays [26] obtained at a lower mass scale. All experimental results are consistent with 
each other assuming the QCD running prediction from RGE. 

The main limitation in the extraction of mb{Mz) from the R^^ measurement is the- 
oretical. If a calculation with resummed LL logarithms [27, 28] could be used, a larger 
range of Ucut could be exploited. This could potentially lead to a lower uncertainty. 

Improvements to the precision of mb{Mz) arc not expected from combining the differ- 
ent measurements because they are largely limited by common systematic uncertainties. 
Other methods will likely be needed at future colliders in order to obtain more precise 
determinations of the 6-quark mass at high energy. This will be important to interpret 
the precise measurements at the Linear Collider in searches for new physics. As an ex- 
ample, a future linear collider operating at ^/s = 500 GeV will produce Higgs bosons 
copiously (if they exist). Since the decay branching fraction into 6-quarks is expected to 
be proportional to the mass squared, measurements of this decay channel would be very 
sensitive to the exact value of the mass at that scale. 
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Ucut 


i?2 


i?2 




i?3 




0.003 


0.505 ±0.003 


0.481 ±0.003 


1.062 ±0.004 


0.342 ±0.002 


0.355 ±0.003 


0.004 


0.553 ±0.003 


0.527 ±0.003 


1.060 ±0.003 


0.329 ±0.002 


0.345 ±0.003 


0.005 


0.585 ±0.003 


0.562 ±0.003 


1.053 ±0.003 


0.317 ±0.002 


0.333 ±0.003 


0.006 


0.611 ± 0.003 


0.591 ± 0.003 


1.046 ±0.003 


0.306 ± 0.003 


0.321 ± 0.003 


0.007 


0.633 ± 0.003 


0.615 ± 0.003 


1.040 ±0.003 


0.295 ± 0.003 


0.310 ± 0.003 


0.008 


0.651 ±0.003 


0.635 ±0.003 


1.040 ±0.003 


0.286 ±0.002 


0.300 ±0.003 


0.009 


0.667 ±0.003 


0.653 ±0.003 


1.031 ±0.003 


0.277 ±0.002 


0.289 ±0.003 


0.010 


0.681 ±0.003 


0.669 ±0.004 


1.027 ±0.002 


0.268 ±0.002 


0.280 ±0.003 



Table 6: Summary of experimental two and three-jet rates, with their total uncertainty, 
as a function of ycut [7] . 



Ucut 


i?4 


i?4 


T3bl 


0.003 


0.1148 ±0.0013 


0.1248 ±0.0009 


0.920 ±0.013 


0.004 


0.0911 ±0.0012 


0.1018 ±0.0007 


0.895 ± 0.015 


0.005 


0.0757 ±0.0012 


0.0856 ±0.0007 


0.885 ±0.016 


0.006 


0.0642 ±0.0011 


0.0729 ±0.0007 


0.882 ±0.017 


0.007 


0.0555 ±0.0010 


0.0628 ±0.0006 


0.884 ±0.019 


0.008 


0.0486 ± 0.0010 


0.0551 ± 0.0006 


0.88 ±0.03 


0.009 


0.0432 ± 0.0010 


0.0486 ± 0.0005 


0.89 ±0.02 


0.010 


0.0380 ±0.0010 


0.0431 ± 0.0005 


0.88 ±0.03 



Table 7: Summary of experimental four-jet rates, with their total uncertainty, as a func- 
tion of ycut [7]. 
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Figure 1: Four-jet rate and its uncertainty as a function of (a) the 6-purity and (b) the 
hght-purity {£ = udsc, Ucut = 0.0065). Chosen working points are marked with arrows, 
and correspond to efficiencies of = 35% and = 86%, respectively. The statistical 
(data and simulation) and total uncertainties are shown. 



17 




" A ^ ^ ^ ^ 
0.005 0.010 .. 



.1.04 



DELPHI 




b) 



0.96 




0.005 



0.010 




0.005 



0.010 



0.005 



0.010 



Figure 2: Comparison between the measured b and i = uds jet-rates and predictions 
from the Pythia 6.156, Herwig 6.2 and Ariadne 4.08 generators, b-d) Ratio of 
data to the different generators. The shaded area shows the one standard deviation 
relative uncertainty (statistical and systematic added in quadrature) of the experimental 
measurement. 
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Figure 3: Comparison between the event-tag (empty circles) and double-tag (full squares) 
techniques for the measured (a) and (b) R^^ observables. The event-tag result of R^^ is 
taken from [6]. The combined statistical (inner bars) and total uncertainty of the exper- 
imental data are shown. The results are compared to the predictions from the Herwig 
6.2 (sohd), Pythia 6.156 (dashed), and Ariadne 4.08 (dotted) event generators. The 
lower insets of the plots show the ratio of data to the different generators. Also shown 
as the shaded area is the one standard deviation relative uncertainty (statistical and 
systematic added in quadrature) of the data. 
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Figure 4: (a) Comparison between the B!'^ measured with a double-tag technique and pre- 
dictions from the Herwig 6.2 (sohd), Pythia 6.156 (dashed), and Ariadne 4.08 (dot- 
ted) event generators. The combined statistical (inner bars) and total uncertainty of the 
experimental data are shown. The lower inset shows the relative deviation of the mod- 
els to the data. Also shown as the shaded area is the total one standard deviation 
relative uncertainty (statistical and systematic added in quadrature) of the data. Be- 
low Ucut = 0.002, the flavour tagging procedure fails and data results from the 1994 
and 1995 data samples are not consistent with each other, (b) Comparison between 
the measured R^^ and theoretical predictions: massive LO predictions and approximate 
(massless) NLO corrections for the pole and running 6-quark mass definitions. Refer- 
ence 6-quark masses were obtained by evolving the average of low energy measurements 
f^biiT^b) = 4.20 ± 0.07 GeV/c^ [14] to the Mz scale as explained in Section 4.5. Hadro- 
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Figure 5: (a) Massive LO results extracted from B!f (data points) compared with the re- 
sult obtained at LO in the analysis [6]: Mb = mb{Mz) = 3.29±0.34 GeV/c^. In the LO 
result, no theoretical uncertainties are shown. Results obtained from R'f using massless 
NLO corrections include theoretical uncertainties estimated as explained in Section 4.5. 
They are shown for the (b) running and (c) pole mass definitions and are compared with 
the results obtained at NLO in the R^^ analysis [6]: mb{Mz) = 2.85 ± 0.32GeV/c2 and 
Mf, = 4.47 ± 0.85 GeV/c^, respectively, shown as ±la shaded band with its central value 
as a dotted line. Predicted values from the QCD calculations at low energy, described in 
Section 4.5, are also shown as solid lines. 
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b-mass measurements at e+e- 




Figure 6: The energy evolution of the MS'-running 6-quark mass mh{Q) as mea- 
sured at Lep. Delphi results from R^^ [6] at the Mz scale and from semileptonic 
-B-decays [26] at low energy are shown together with results from other experiments 
(Aleph [3], Opal [4] and Sld [5]). The masses extracted from LO and approximate 
NLO calculations of i?4^ are found to be consistent with previous experimental results 
and with the reference value mb{Q) (grey band) obtained from evolving the average 
TTibimb) = 4.20 ± 0.07 GeV/c^ from [14] using QCD RGE (with a strong coupling con- 
stant value a,{Mz) = 0.1202 ± 0.0050 [25]). 



